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Abstract. - We introduce a finite-connectivity ferromagnetic model with a three-spin inter- 
action which has a crystalline (ferromagnetic) phase as well as a glass phase. The model is not 
frustrated, it has a ferromagnetic equilibrium phase at low temperature which is not reached 
dynamically in a quench from the high-temperature phase. Instead it shows a glass transition 
which can be studied in detail by a one step replica-symmetry broken calculation. This spin 
model exhibits the main properties of the structural glass transition at a solvable mean-field 
level. 



Over the last decades there has been a growing interest in understanding the structural 
glass transition in complex materials. A major theoretical role has been played by mean-field 
theories based on fully connected p-spin glass models [0-1). In such models, it was possible 
to identify a purely dynamical transition which accounts for the off-equilibrium behaviour of 
the models, and which is believed to present some analogies to the structural glass transition 
in realistic systems 

Although this analogy is very fruitful, one must keep in mind that the mean-field spin glass 
models which have been considered are rather remote from the structural glasses. Let us men- 
tion just the three most obvious differences: they have no crystalline state, their Hamiltonian 
has quenched-in disorder, in the form of an extensive number of quenched coupling constants, 
and these couplings are of infinite range. A lot of work has already been devoted to cure these 
defects, in order to get a much closer analogy between some lattice spin models and structural 
glasses. This should allow to clarify some basic issues as for instance the presence of hetero- 
geneities or the role of geometrical frustration. A few years ago, it has been realized that the 
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existence of quenched disorder is not a necessary ingredient for a system to have a spin glass 
phase: several mean-field spin models were constructed which display a discontinuous spin 
glass transition typical of the p-spin glasses |^ . More recently, numerical simulations of some 
unfrustrated 3-dimensional spin problems, with purely ferromagnetic multi-spin interactions, 
have shown a glassy behaviour which persists for very long times 1^ . A similar behaviour was 
found also analytically in 2 dimensions |lO[ . 

In this paper, we introduce a similar ferromagnetic model having multi-spin interactions 
and finite connectivity. The model is defined on a Bethe hyper-lattice (or Husimi tree, the 
analog of a Bethe lattice for systems with plaquette interactions), where it can be solved 
completely using the replica/cavity method. This model thus provides the correct mean-field 
(Bethe-) approximation for these ferromagnetic multi-spin systems. As in the usual Bethe lat- 
tice, it is a finite-connectivity system (each spin interacts with a finite number of neighbours), 
but has a locally tree-like structure which makes it solvable at the mean-field level. The spin- 
glass version of the model, with random boundary conditions, was considered first in Our 
solution shows a phase diagram qualitatively similar to that of usual p-spin models, with a 
static and a dynamic glass transitions appearing at two different temperatures. On top of the 
glass phase, it also displays a purely ferromagnetic phase which is the analog of the crystalline 
phase in supercooled liquids. 

The model is defined by the Hamiltonian 



where Si — ±1 (i = 1, A^) are Ising spins. The model is ferromagnetic: we take Jijk ~ 1 
for all the triples k] belonging to the set E of hyper-edges (plaquettes) . So the system is 
not frustrated and the configuration Si = 1 (i = 1, .., A^) is a ground state for every possible 
choice of E and satisfies simultaneously all interactions. The plaquettes are chosen randomly, 
according to the two following ensembles: 

• In the first ensemble, the only constraint is that each spin belongs exactly to fc -I- 1 
plaquettes. The set of plaquettes builds up a hyper-graph which is locally a Husimi 
tree of fixed connectivity A; -I- 1. This structure of hyper-edges eventually loops back 
onto itself, but the typical length of the loops scales like ln(7V). This hyper-graph was 
used analogously in frustrated spin models studied in the context of spin glasses p^ . 
In particular one should notice that its structure is not disordered on any finite length 
scale: disorder comes in only through the loops which have diverging length scale and 
so the system is only very weakly disordered. Basically this construction amounts to 
having a Cayley hyper-tree with random hyper-edges closing the boundary. 

• Our results can be easily generalized to hyper-trees with fluctuating connectivities. Such 
models, with Poisson-distributed connectivities, were introduced recently in ||l3|,|lj] as 
a simple case of an optimization problem. As in the case of ferromagnetic systems with 
multi-spin interactions in finite dimensions |l0| , it was observed that purely ferro- 
magnetic interactions may lead to an effective dynamical frustration. This was shown 
in Jl^ by comparing lowest metastable states in a ferromagnetic finite-connectivity 3- 
spin model with the ground states of the corresponding frustrated spin-glass version 
{Jijk = ±1 randomly), which were found to share exactly the same statistical features. 
As also seen in Q by analyzing the off-equilibrium low temperature dynamics of both 
models, the ferromagnetic system is not able to equilibrate due to the entropic domi- 
nance of metastable states. 
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The results can also be generalized to interactions between clusters oip spins (p > 3), also with 
spins belonging to clusters but interacting pairwise, which can be viewed as a representation 
of geometrical frustration. 

At first we concentrate on an analytical approach to the equilibrium behaviour based on 
the replica trick The results will be compared later on with data obtained form Monte 
Carlo simulations. The free-energy density can be read off from the n ^ limit of the n-fold 
replicated free energy at temperature T = f3~^ . The formulas are written here for a general 
hyper-tree with fc + 1 hyper-edges per site, we shall then discuss in more details the case 
fc = 3. This is the lowest value of fc for which ferro-magnetic and glassy behavior are present, 
for k — 2 the system is paramagnetic at all temperatures pT| . 

The natural order parameter in the replica approach is a probability distribution c(it) 
on the set t? £ {±1}", which counts the fraction of sites i having replicated spin Sj^ — (j^ 
(a = 1, n). For details on this approach sec Rcfs. In terms of this order parameter, 

the replicated free energy reads 



-/?/„ = extr 



k+1 



In [ c(ai)'^c(a2)'=c(a3)'e'5S"=i<"=°'^3 j _ fcln ( ^c(a)'=+M 



(2) 



We notice that if wc had considered spin-glass couplings Jijk = ±1 with symmetric probability, 
we would have found cosh(/3^"^-^ cr^'crfCTf ) instead of the exponential in (Q). 

One needs to find a distribution c{<7) which makes this free energy stationary, which means 
that it must satisfy the saddle point equation 

c(a) cx ^ c(a2)^c(a3)' exp [/sf^ aV^VgM (3) 

(72, (73 \ a—1 / 

It turns out that there exist three solutions, corresponding to a paramagnetic, a ferromagnetic 
and a spin-glass phase of the system. 

The simplest solution is the paramagnetic one, Cpm{S) = 1/2" for all a. This solution is 
expected to be valid at sufficiently high temperature. The paramagnetic free-energy density 
is given by 

-/3/pm = In 2 + ^ In cosh(/3) . (4) 

At lower temperature there appears a ferromagnetic solution. There, replica symmetry 
(RS) is expected to hold, and c(a) depends on a only via X^a"'"- Furthermore there are no 
site to site fluctuation since all nodes of the hyper-tree are equivalent, so one can look for a 
ferromagnetic solution of the type: 

= (2cosh/?fe)» • 

This is indeed a stationary solution provided the effective magnetic field /i is a solution of the 
equation: 



tanh(/3/i) = tanh(/3) tanh(/3fc/i)^ 



(6) 
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Nontrivial solutions (with h ^ 0) exist only for temperatures below a certain Tms, and disap- 
pear abruptly above Tms- The resulting free energy reads: 

-Pffrn = ln2 + In [cosh(/3) cosh{pkhf + sinh(/3) smh{l3kh)^] - fcln [cosh(/3(fc + l)/i)] 
3 

(7) 

The ferromagnetic transition is of first order: one needs to find the largest, i.e. locally stable, 
of the two nontrivial solutions for h. For temperatures slightly below T^s the free energy 
of this ferromagnetic solution is larger than that of the paramagnetic one and the latter 
stays globally stable. Only at a temperature Tj™ < Tms, the ferromagnetic solution becomes 
thermodynamically dominant (lower free energy), inducing a first order transition. For the 
case fc = 3, one gets Tms — 1-63 and Tfm — 1.21. 

A third solution of (||) corresponds to a glass phase. As we shall see, this phase is easily 
accessed by Monte Carlo simulations with decreasing temperature. Its physical origin lies in 
the fact that if one of the spins on a plaquette is down, then the effective interaction among 
the other two spins in the same plaquette becomes antiferromagnetic. Therefore there is no 
restoring force towards the ferromagnetically ordered state. In order to find the glass phase 
in the replica approach, one has to look for a c{a) which breaks replica symmetry. 

We have found a one step replica-symmetry breaking (RSB) solution, which is even with 
respect to separate changes of sign of all its variables. This solution is somewhat involved 
because we are dealing with a finite connectivity system. We need to build up, on each site 
i, a distribution of local fields Pi{h). The local fields on site i in the various pure states a 
are iid variables chosen from the distribution Pi{h). This distribution fluctuates from site to 
site and the order parameter is a functional: the distribution of the functions Pi. In replica 
language, this one-step RSB amounts to a solution Crsb{'^) of the saddle point equations with 
the following structure: 



/n/m . 
dA^(A) J]^ / dMa'^(ua|A)exp /3ua ^ o- 
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6=(a-l)m+l 



(8) 



It depends on the real (for n ^ 0) number m, on the function /i which is a probability 
distribution, and on the functions (/)(m|A)(2 cosh(/3u))™ which are probability distributions 
on u conditioned to a given value of A. The link between the replica approach and the 
distributions Pi (K) is reviewed in |l^ . 

Despite difficulties in finding the saddle point, it is easy to see that the spin-glass solution 
is stable against ferromagnetic fluctuations. Indeed, the explicit computation of the replica 
Hessian matrix corresponding to free energy (^ shows that if cia) is even, then the matrix is 
positive definite on the subspace of functions via) which are odd in at least one variable |(^ . 

We have used the method recently proposed in ||l^ which allows to determine the sad- 
dle point by a population dynamics of local fields. We have found that a nontrivial solu- 
tion appears for m = 1 below a certain dynamical temperature Tc. This solution appears 
discontinuously (the fields are not small close to the transition). Below a temperature Tk 
this solution becomes thermodynamically relevant (letting aside the ferromagnetic state) and 
the saddle point value of m is smaller than one (see inset of Fig. |l|). Between Tk and Tc 
the parameter m sticks to one and the equilibrium free energy is that of the paramagnetic 



(i)The result holds more generally for diluted p-spin models for p > 3, while for p = 2, the spin-glass solution 
exists, but is unstable against ferromagnetic fluctuations. 



S. Franz , M. Mezard , F. Ricci-Tersenghi , M. Weigt and R. Zecchina 



5 



phase. This intermediate phase corresponds, as usuaUy, to broken ergodicity with exponen- 
tially many ergodic components, i.e. extensive configurational entropy. For the case A: = 3 we 
have found Tc = 0.745(5) and Tk = 0.660(5), while the configurational entropy at equals 
Sconf{Tc) — 0.063(5) and, as it should, vanishes at Tk- 

A simple analytic approximation to this one-step RSB result can be obtained by a varia- 
tional approximation. In replica language, one can use for instance an approximate form for 
the order parameter c{a) of the type Q 



^ ' -I J- J Du{2 cosh 13 Au)'^ ^ ' 



a=l 



where Du = du exp{—u^ /2) / \/2tt . The free-energy density resulting from this Ansatz has to 
be optimized with respect to the variational parameters A and to. Wc find T^"'*^ — 0.752 
and T]^'^'' — 0.654, which coincides well with the numerical solution of the exact saddle-point 
equation. 

The temperature Tc is a dynamical temperature where the relaxation times diverges and 
the system becomes non ergodic, while the temperature Tk is the Kauzmann temperature 
where the system has a thermodynamic phase transition (as in many other cases the fact 
that Tc is different from Tk is made possible by the mean field nature of the problem 0). 
This spin glass solution presents all the properties of the discontinuous spin glasses in which 
one-step RSB is exact [Q. However, we can not exclude a second static transition with full 
replica symmetry breaking at a temperature even smaller then Tk, as it happens in the fully 
connected Ising p-spin model |T7[ . 

Note also, that the RSB solution for c{(j) is even under reversal of all n spins, and, thus, 
is also valid for the spin-glass version of the model, where couplings are set to ±1 randomly. 
The glassy and the paramagnetic phase turn out to be indistinguishable in both models. 

In order to clarify the physical implications resulting from the existence of three saddle 
point solutions, we have performed Monte Carlo (MC) simulations. First of all we have 
determined Tms as the maximum temperature at which the ferromagnetic state is locally 
stable. In Fig. |l| the heating curves show the results of MC simulations ~ started a T = 
with all spins up - during which the temperature is slowly increased. The static transition 
temperature Tfm is determined as the point where the time spent by the simulation in the 
paramagnetic and in the ferromagnetic state are equal 

If we start, however, from a random configuration, the evolution of large systems {N = 
99999 in our simulations) is completely insensitive to the presence of a ferromagnetic ground 
state, the time to find it being exponentially large in the system size |T4| . The system remains 
unmagnetized even if quenched below Tj:„n and may thus undergo only the glass transition. 
We have measured the stationary spin-spin correlation function, and have determined the 
critical point Tc from the divergence of the correlation time. The result Tc = 0.75(1) is 
perfectly compatible with our analytical estimates. In this respect, the 3-spin ferromagnetic 
model with finite connectivity is completely different from any 2-spin unfrustrated model or 
any p-spin unfrustrated fully connected one. Indeed, in the latter cases the stability argument 
presented above does not hold and a simple coarsening dynamics pushes the system towards 
the unfrustrated ground state |p)| , thus ruling out the glassy off-equilibrium behavior discussed 
above or in 

In Fig. nTwe show the results of a set of cooling experiments on the unfrustrated model 
(identical curves have been obtained for the frustrated model). When the simulation goes 



(^)The same trivial dynamics applies to any unfrustrated model whose connectivity diverges with A'^, because 
thermal fluctuations in the magnetization are amplified, until the ground state is reached. 
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Fig. 1 Fig. 2 

Fig. 1 - Average energy as a function of the temperature for the ferromagnetic model with fe = 3. 
The rightmost vertical lines are the analytic prediction for the ferromagnetic transitions, while the 
leftmost vertical arrows are those for the spin glass transitions. Continuous lines: results of simulated 
annealing experiments (cooling rates from 10^ to 10^ MCS per AT — 0.01), where the magnetization 
stay zero and the energy converges to the spin glass threshold energy. Squares: static energy in 
the glassy phase. Crosses: the system is initialized fully magnetized and then slowly heated. Data 
are indistinguishable from the superimposed analytic curve. INSET: saddle point value for m as a 
function of temperature. The line is the fit 1.41T + 0.138T^ to T < Tk data. 

Fig. 2 - Continuous (dashed) lines represent the static (dynamic) transition lines for a 3-spin model 
with average connectivity 37 calculated with a variational Ansatz. Upper (lower) lines refer to the 
ferromagnetic (spin glass) transition. Crosses (with errors): estimations of the critical lines from 
Monte Carlo simulations. Squares: results of the algorithm which gives the exact 1-RSB solution. 
Black dots on the T = axis mark the exact results for the ferromagnetic model. 



through the ferromagnetic critical points (the 2 rightmost vertical lines) the system's evolution 
is completely unaffected, and the average magnetization stays near zero. The relaxation 
process is in fact strongly slowed down only when the spin glass critical points are reached 
(marked by the 2 leftmost vertical arrows). Note also that below the spin glass critical point 
the energy relaxation is almost absent and the asymptotic energy strongly depends on the 
cooling rate. Indeed the most effective relaxation is the one happening close to the critical 
point. All the above features are typical for structural glasses. 

We have checked that some amount of disorder does not destroy the above picture, by 
studying a system with fluctuating connectivity, in which the set E of hyper-edges contains 
7A'^ randomly chosen triples fc], i < j < k. In this case the average connectivity c = 87 
can be varied continuously and critical lines in the (7,r) plane can be evaluated. Fig. ^ 
summarizes our analytical and numerical findings on this model, which are in full qualitative 
agreement with those obtained in the case of a hyper-tree with fixed connectivity. 

The analogy with supercooled liquids can be pushed forward once we identify the fer- 
romagnetic ground state with the crystalline phase, the paramagnetic state with the liquid 
phase and the spin glass state with the glassy phase. If we do not force the system towards 
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the ferromagnetic ground state, below Tfm the system remains in the paramagnetic phase 
and it is thus supercooled. Decreasing further the temperature the system undergoes a glass 
transition, even if no quenched frustration is present. The frustration is self-induced by the 
dynamics: being unable to find the ferromagnetic ground state, the system is typically in 
a configuration where a finite fraction of the interactions are unsatisfied and no long range 
ferromagnetic order arises. Indeed, the ground state can not be found (in polynomial time) 
by simply exploiting local information, the only one available when a spin is updated in a 
single spin flip dynamics. We expect a long lived glassy regime to be present also in short 
range versions of our model, where it should be destabilized by slow enucleation processes, as 
already seen in numerical simulation of some three dimensional systems |^ . 

Spin glasses were originally understood as intrinsically disordered and frustrated systems. 
A few years ago it was found that the disorder is not really necessary to induce a spin glass 
phase y . We have shown here, studying diluted mean-field models, that quenched frustration 
is not necessary either. Our findings give a strong support to the use of disordered and 
frustrated models (e.g. p-spin glasses) to describe structural glasses, which are by constitution 
neither disordered nor frustrated. 

* * * 

M.W. acknowledges the hospitality of the ICTP and S.F. that of the LPTMS. This research 
has been supported in part by the SPHINX project of the European Science Foundation. 
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